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I. INTRODUCTION 



Inflationary paradigm besides remarkable successes in describing the cosmological data 
and CMB temperature anisotropies [sl , has a very appeahng theoretical and model building 
feature: inflation considerably relaxes the dependence of late time physics on the pre- inflation 
initial conditions ji]. A particular set of initial pre-inflation conditions which will be of our 
interest in this work is the anisotropic but homogeneous initial conditions for metric, the 
Bianchi type cosmology jsl]. 

Although not directly related to inflationary models, Wald proved his elegant cosmic no 
hair theorem In the presence of a cosmological constant, anisotropic but homogenous 
deviations from de Sitter inflating background with energy momentum tensor respecting 
strong and dominant energy conditions will be exponentially damped by the dynamics of 
the system within a few Hubble times scale. Wald's cosmic no hair theorem has been the 
basis of the arguments on how inflation washes away the anisotropies, regardless of the 
details of the inflationary model in question. This expectation is named the cosmic no-hair 
conjecture. 

Inflationary models are generically using scalar flelds as inflaton(s) and, dealing with 
scalars it is straightforward to check that cosmic no hair conjecture is valid in these models. 
For an early conflrmation of this conjecture in the context of chaotic inflationary models see 
7|.) This may be seen from the fact that the amplitude of vector perturbations in standard 
cosmic perturbation theory are (exponentially) damped by the inflationary expansion of the 
background ^. 

Presence of vector flelds, however, may change this picture. This may happen, for exam- 
ple, in the vector inflation models, where there are vector flelds turned on in the background 
a. inflaton fl Alternative., o.e „,ay add U(l) gauge field, with a .pec.fical.v ta,- 
lored kinetic term to a standard scalar driven inflation model. The kinetic term can be 
appropriately chosen such that it can compensate the exponential damping of the vector 
modes HI]. In this latter class of models the modifled Maxwell kinetic term provides the 
setting for violation of assumptions of Wald's cosmic no hair theorem, allowing for growth 
of the anisotropies. Thus, there exists a counter example to the cosmic no-hair conjecture. 
Nonetheless, demanding a successful inflationary model restricts the amplitude of anisotropic 
^ We note that vector inflation models, due to broken gauge invariance, generically have ghost problem, 
inducing instability in the background inflationary dynamics 2, lo| . 
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Bturbations to remain small Il2| . compatible with the current observational bounds 
. (Anisotropic inflationary models have also been discussed in 14j-|l6|.) 
In two of us introduced a novel inflationary scenario, non-abelian gauge field inflation 
or gauge-flation for short. In this model, which will be briefly reviewed in sectionllTl inflation 
is driven by non-Abelian gauge fleld minimally coupled to gravity^ It was shown that 
non-Abelian gauge fleld theory can provide the setting for constructing an isotropic and 
homogeneous inflationary background |2|. Despite of turning on background gauge flelds, 
the isotropy of the background was reinstated by using the global part of the gauge symmetry 
of the problem and identifying the SU{2) subgroup of that with the rotation group. We 
argued that this can be done for any non-Abelian gauge group, as any such group has an 
SU (2) subgroup. Therefore, our discussions can open a new venue for building inflationary 
models, closer to particle physics high energy models, where non-Abelian gauge theories 
have a ubiquitous appearance. 

Due to the gauge-vector nature of our inflaton fleld a question that may arise naturally 
is the classical stability of gauge-flation against the initial anisotropies. In fact, the rotation 
symmetry of the universe was essential to set a consistent ansatz with an isotropic universe. 
Hence, it is not obvious if the special conflguration adopted in the previous work is stable 
against nonlinear perturbations in the anisotropic background universe. (Stability against 
linear perturbations around FLRW gauge-flation background was shown in ^ in the context 
of cosmic perturbation theory.) Moreover, one can ask if gauge-flation model is conformed 
to the cosmic-no-hair conjecture, and if speciflc observational features related to anisotropy 
originated from the gauge flelds can be observed. In this paper, we will tackle these questions. 
Explicitly, we show that starting with generic anisotropic initial conditions within Bianchi 
type-I cosmology, dynamics of gauge-flation suppresses the anisotropies. That is, isotropic 
FLRW cosmology is an attractor of the gauge-flation dynamicsjf] Since the anisotropies 
are damped and basically washed away in the flrst few e-folds of the inflationary period, 
gauge-flation predicts a very small (negligible) statistical anisotropy in the power spectrum 
of CMB fluctuations. 



^ The gauge-flation respects the gauge and difl'eomorphism invariance and hence is free of the ghost problem 
of vector inflation models. 

Presence and effects of non-Abelian gauge fields in the context of inflationary models have also been 



studied in 



14| and [17|. These are different than the gauge-flation case because there is also a scalar 



(inflaton) field in these models. Moreover, the model discussed [14] involves an action with SU{2) local 



in |l7| h 
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gauge symmetry, while the one in |17j involves an action with SU (2) global symmetry. 



This paper is organized as follows. In section |TT1 we introduce the gauge-flation setup, 
and briefly review the model for isotropic FLRW cosmology and inflation in this model. 
In section IIIH we consider a Bianchi type-I background metric and, using quasi-de Sitter 
approximations we study the background anisotropic inflation in gauge-flation setup. Prob- 
ing space of solutions and classical trajectories of the model, we show that regardless of 
the initial value of anisotropics, our initially anisotropic inflation evolves toward isotropic 
solution within a few e-folds and after that it effectively mimics the behavior of the isotropic 
quasi-de Sitter inflation. That is, the isotropic background is the attractor of the system 
with generic anisotropic initial conditions. We backup the analytical results of this section 
by numerical analysis over a range of initial-values. Finally, in section llVt we summarize 
our results and make concluding remarks. 



II. GAUGE-FLATION SETUP 



Gauge-flation is based on a non-Abelian gauge theory minimally coupled to Einstein 
gravity. In [l, 2| it was shown that Yang-Mills action for gauge theory part cannot lead to 
accelerated expansion and one should consider addition of other terms. One particularly 
convenient choice of the action is 



S 



d X\ 



-9 
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(II.l) 



where we have set SvrG = Mpj^ = 1 and is the totally antisymmetric tensor. The gauge 



field strength F"'^^ is given by 



p^a 



A'^ ^ A°- — Of:"- A^ 



(11.2) 



where A"-^, /i = 0, 1, 2, 3 label the spacetime indices and a the internal gauge indices. We 
choose the gauge group to be SU{2) and hence a = 1, 2, 3. 

Since yf—ge^^^'' is independent of the metric, the metric dependence of the term 
appears only through deig. Its contribution to the energy momentum tensor, denoted by 
T^^, is of the form 

T;y = -p,{t)g,y. (II.3) 

In the absence of the Yang-Mills term, equation of motion for the gauge field gives a constant 
Pk and hence, if k > 0, we have a de Sitter expansion phase. Addition of the Yang-Mills term, 



however, changes this behavior and eventually ends inflation. (Note that energy momentum 
of the Yang-Mills term T^J^ is traceless and hence cannot sustain accelerated expansion 
phase.) 

To be more precise, defining the isotropic expansion rate H as 

3« 



where are the space-like constant comoving time hypersurfaces, we can define the param- 
eter 

which is nothing but the slow-roll parameter in the standard slow-roll inflationary scenario. 
This parameter describes to what extent the universe is close to de Sitter spacetime. In this 
paper we will use the expression "quasi-de Sitter" inflation for denoting a background with 
e ^ 1. Therefore, quasi-de Sitter captures the standard isotropic FLRW slow-roll inflation 
as well as allowing for anisotropic cases. 

From the above argument, demanding quasi-de Sitter inflation (e ^ 1) is equivalent to the 
condition ^y*^ ^ 1 (where pyu is T^y^ projected on the vector normal to Sf). Namely, when 
the contributions to the energy density dominate over the Yang-Mills contributions 
PyJ^^ , the universe is quasi-de Sitter. The time evolution of the system then slowly increases 
pyj^j with respect to p^, and when p^^^ ~ p^, the quasi-de Sitter inflation ends. 

The above analysis can be made more explicit if we consider an isotropic homogenous 



Here, we review some of the 



background. This was the case studied in some detail 

ca.cu.at,o,. of 3 as a wa^up and also fo. fixing scne convention, and notation, we will 

use in the anisotropic case of the next section. We start with the background FLRW metric 

ds^ = -dt^ + a^{t)5ijdx'dx^ . (II. 5) 

Choosing the temporal gauge A"'^ = 0, we can set the ansatz 

A% = = amm, (11.6) 

where are the triads on constant t slices, is consistent with the dynamics and identifies 
the combination of the gauge fields for which the rotation symmetry violation caused by 
turning on vector (gauge) fields in the background is compensated (or undone) by the gauge 
transformations, leaving us with a rotationally invariant background. Note that under both 



of 3d diffeomorphisms and gauge transformations iplt) behaves as a genuine scaler field, 
nonetheless it appears that equations take a simpler form once expressed in terms of (j) 

= ait)m , (11.7) 

or equivalently A"-- = (f){t)6°'^. 

Evaluating f lll.ip for the FLRW metric and the ansatz flll.6p . we obtain 

r 3,0^ g'<P\ ^?W, ,,,,, 

/:.e.= 2(^-^ + -^^), (n.8) 

which is the Lagrangian governing dynamics of the system in the homogenous-isotropic 
sector 21]. Using Lred one can compute the energy density p and pressure P 

P = Pym+ PK, P = ^PyM - PK , (11.9) 

where pym and are respectively contributions of Yang-Mills and k terms of the action 



Note that in order to have positive energy we consider positive n values, hence Py^^ and p^ 
are both positive quantitiesc 

The Einstein equations, the Friedmann equation and the evolution equation for H, are 
then obtained as 



^ = 2^^ + ^ + ^ = + ^) • ("-^'^ 

Using the Friedmann equations ( 111. lip and definitions ( lll.iop and ( 111.4^ . we have 



2Py„ ^jj_^2) 



Pym + P« 

Here, we can see that in order to have quasi-de Sitter inflation the K-term contribution p^ 
should dominate over the Yang-Mills contributions Pym^ i-^- Pk ^ Pym- 
The Yang-Mills equations reduces to a single equation for0 field: 

1 + K^-^)^ + (1 + «:^)^ + (1 - 3«:^)^ = . (11.13) 



^ Note also that the non-Abehan character of the fields, besides in resolution of anisotropy of the background 
caused by turning on the gauge fields, also shows up in the fact that the K-term gives a nonvanishing 
contribution to the energy density p, which is essential for having inflation. 
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In the K dominant limit, i.e. when e <^ 1, the above equation becomes 



2, 




-0, (11.14) 

where ~ means equahty to the leading order in the parameter e. This implies ~ constant 
and then H ~ constant. In this limit, by integrating Eq.( ]II.14|) . we also have a relation 
(j) oc exp Ht oc a. Hence, ip must be almost constant during quasi-de Sitter inflation. Since 
we also have a relation ~ we obtain the formula 

e~V^2(7 + l), (11.15) 

where 7 is defined as 

7 = ^- (11.16) 

Thus, to have a sufficient e-folding number of inflation, we need sub-Planckian field values, 
ip '^l. (Note that we are working in Mpi = 1 units.) 

To summarize, assuming ^ PvAn have a successful isotropic quasi-de Sitter inflation 
period, during which H and ip are almost constant. Since e oc ip'^, the quasi-de Sitter 
condition demands sub-Planckian field values. The time evolution will then increase PyJ^^ 
with respect to p^, and when p^^^ ~ p^ the slow-roll inflation ends. Noting that p -|- 3P = 
'^(Pym ~ Ph)^ inflation (accelerated expansion phase) ends when p^^,^ > p^. 



III. GAUGE-FATION AND COSMIC NO-HAIR 



In this section we study gauge-flation in a homogenous but anisotropic background. In 
this way we examine the generality of the isotropic FLRW gauge-flation. Here, for practical 
reasons we consider gauge-flation in an axially symmetric Bianchi type-I setup but we believe 
that our results are valid for more general anisotropic cases. 

Bianchi type-I axially symmetric metric is described by the line element 

ds^ = -N^dt" + e^^W (e-^-^Wrfx^ + e^'^^'\dy^ + dz^)), (III.l) 

where N is the lapse function, e""^*-* represents the anisotropy and, e"*^*-* is the isotropic scale 
factor. Given the symmetries of the metric, as before we choose the temporal gauge for the 
gauge fields = 0, and make the following modification to the isotropic ansatz (III.6P 

/l^ = e>„ ^1 = ^, ^2 = ^3 = AV^, (III.2) 
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where ibi act as two scalar fields and e° are the triads associated with spatial metric. The 

n 

gauge field one- form of the above ansatz is hence of the form|j 

. . . e-.-«mr... . e-.-...A,<).(«)(r^# + T'..) . ,111.) 

where T° are the SU{2) gauge group generators, [Ta,Tb] = ie'^^^^Tc. It turns out that, similar 
to the isotropic case (111.7p . the equations take a simpler form once written in terms of 

cPit) = a{t)m , (III.4) 

where 

a{t) = e" , (III.5) 

is the isotropic scale factor. 

Inserting the axi-symmetric Bianchi metric and the gauge field ansatz into the gauge- 
fiation action (III. II) . the total reduced (effective) action is obtained as 

a^r „.o .o.„ (2 + A6)02 , (A-4(A6 - 1)02) (i + 2A6)02 



Sred — N 



-3a' + a^3+ '-^^) + a ' \^ ' + 



\^ a^' a2 2A4 



(A6 - 1) A 00 , (2 + A6) A2 02 (2 + A^) g'cj^' , 3 ng'cj^^ cj^^' 



+ 2^ ^-im + ^ — \ 



(III.6) 



A4 Aa2 A4 A2a2 2A2 2 _ 

where a dot denotes derivative with respect to the time coordinate t. As we see, the above 
action depends only on a and not a. Therefore, momentum conjugate to a is a constant of 
motion. This constant may be chosen on the physical basis that the anisotropy a should 
vanish for the isotropic gauge field, i.e. when A^ = 1. With this initial condition we obtain 

^ = 7 ^ ■ 

2a2(3 + A-4(2 + A6)g) 

We can deduce the field equations corresponding to 0, A and a from the above effective 
action. However, instead, we find it more convenient to determine the energy-momentum 
tensor and study the Einstein equations, as well as the field equation for A. Substituting 
the ansatz ( 1111.31) and the metric ( IIII.II) into T^^, for the gauge fields, we obtain a diagonal 



homogenous tensor, 

T; = diag(-p(t),Pi(t),P2(t),P2(t)). 



A similar ansatz has been considered in 



8 



One can decompose the energy density p as 



(III.8) 



where and p^^^ are respectively contributions of F*^ and Yang-Mills terms given by 



Pk „ 4 r, : 

2 a* a"^ 



Pi 



(III.9) 
(III.IO) 



Note that p^ is only a function of and not A. As mentioned before, we consider positive 
K, values, so p^ and p^^ are both positive quantities. Furthermore, Pi and P2 are 



Po = P+ ip, 



where the isotropic pressure P and the anisotropic pressure P are given by 

P=-PK+ ^Pym 



(iii.ii) 

(III.12) 



(111.13) 



A 



o + -)-J -. 



-A^^a ^ ' \' a A^ ^A ' ''-a ^ ' \' a 

As we see, in the isotropic case A^ = 1 (6" = and A = 0), P vanishes and hence Pi = P2 = P. 
The independent gravitational field equations are 



•2 -2 P 

« -a =-, 



ci + 36" 



P2 - Pi P 



3 3 
_ 3p + Pi + 2P2 _ 2 

6 " 3 



o PyM • 



(III. 14) 

(III.15) 
(III.16) 



Note that (1III.15I) implies that the evolution of a depends on the anisotropic part of the 
pressure P, i.e. P is the source for the anisotropy and in the absence P the anisotropy a 
is exponentially damped, with time scale oT^ . 



A. Analysis in Quasi-de Sitter regime 

So far we computed T^,^, the reduced action and a without assuming the K-term dom- 
inance. At this point, we simplify and analyze the equations assuming quasi-de Sitter 
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inflation, in the sense that the parameter e 



a 



is small during inflation. Note that we only impose the quasi-de Sitter condition on the 
isotropic parts of the expansion and the rest of variables are not enforced to have slow 
dynamics. 

From the combination of (1111.14^ and (1111.16^ we obtain e in terms of p^, Pym ^ 



2p.M + 



(III.17) 



P« + Pym + 3(t2 ' 

which implies ^ p^j^^ and p^ ^ for having a very small e. Therefore, from (IIII.14I1 
and after using ( 1111.91) . we obtain 

,2 i4 12 



P. = ^^^ = 3d^ + 0(e). (III.18) 
As we have learned in the previous section, the constancy of implies 

(111.19) 

Moreover, noting that P is a term from the contribution of Yang-Mills part, and the corre- 
sponding isotropic pressure is also positive, we learn that p^ ^ P. Then, using fllll.lSP we 
obtain a& ~ P which gives 

^2 (111.20) 

and as P is a quantity at most of same order as Pym, we obtain &^ ^ PyAf ^ result, 
during quasi-de Sitter inflation, we have the following approximation for e 

e ~ ~ (111.21) 

p^ 

where ~ means equality to the leading order of e. Using the above relation and fllll.l9p . 
ignoring the a terms in p^^^ , we flnd 



'^2 



3a2A4 



(2 + X')iX'l + ^) - ^(1 + A^) + (1 + 2A6) 



(111.22) 



where 7 = fs^- As we will show through analytical calculations, ^ is at most an order 
one quantity. Therefore, in order to have a successful quasi-de Sitter inflation (e ^ 1), we 
should have 

— < 1 , 4i < 1 ' (™-23) 
a aX^ 
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for large and small A, respectively. In other words, recalling ( ]III.2p . similar to the isotropic 
inflation, our field values ipi^s should have (physically reasonable) sub-Planckian values dur- 
ing the quasi-de Sitter inflation. 

From a& ~ P and P < Pym, we obtain 

•2 / r>\^ / \2 



That is, that the value of ^ in our anisotropic gauge-fiation model cannot be more than e^. 
To summarize, so far we have shown that during the quasi-de Sitter inflation still means 
dominance. Precisely, e ~ -^-f^ and which represents the anisotropy of the expansion 

of the universe, can be at most of order e^. To proceed further we need to analyze dynamical 

field equations. 

Using (lIII.14p . (IIII.2ip and ignoring a terms, we obtain the following relation to the 
leading order of e, 

3d' -p., (III.25) 
Neglecting a terms, we can deduce the field equation corresponding to A as 

(2 + A^)(AA + aXX + 2^AA) - eA^) ^ + A'(A^ - 1)(^ + ^ + ^) = 0(111.26) 

which, using (1III.19P and keeping the leading orders, can be simplified to 

(2 + A^)(AA + 3dAA) - 6X^ + X\X^ - 1)(2 + A^)^' ~ 0. (III.27) 

One can see that A = is a singular point of the above equation and that the dynamics 
do not mix A > and A < regions. That is, if A is positive (negative) initially, it always 
remains positive (negative) during the inflation. 

Since (1III.27P is a nonlinear second order differential equation which has no explicit time 
dependence, its solution will be of the form A = A(A) and hence 

A = ^A. (III.28) 
dX 

In terms of derivatives with respect to dN = adt, and denoting ^ by a prime, we obtain 

dX' 1 ,6A' A(A6-1), ,2,, ,,,, , 

which implies that A'(A) is an odd function of A, A'(— A) = — A'(A). 
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Using Mathematica, the above equation can be studied by the phase diagram method, and 
in Figure [H we have presented the behavior of the solutions in the A' — A plane. Apparently, 
all of trajectories approach to the isotropic fixed point A^ = 1. Next, we give the asymptotic 
analysis to confirm that the isotropic inflation is an attractor in the phase space. 




-10 -5 5 10 0.90 0.95 1.00 1.05 1.10 



FIG. 1. The phase diagram in the A' — A plane, the vertical axis is A' and the horizontal is A. Both 
figures show existence of attractor at |A| = I, corresponding to the isotropic FLRW background. 
The left figure shows the phase diagram over a large range of values for A, while the right figure 
shows the phase diagram for A in the vicinity of the attractor solution A = 1. The left figure 
explicitly exhibits the A'(A) = — A'(— A) symmetry. 



B. Asymptotic Analysis 

To study the behavior of the system more precisely, it is convenient to rewrite ( 11112711 as 
(2 + A*^)(^ + 3^) - 6^ + (A*^ - 1)(2 + A^) ^ , (11130) 
in which we distinguish three different regions for the value of A: 
I) A^ — 1 close to zero (A^ ~ 1), 
II) A^ in the vicinity of zero, 
III) large A^ limit. 
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Furthermore, defining 



E = - = a', (III.31) 

a 



and using (IIII.19I) in flIII.7|) . we can write S as 



We will solve ( 1111.301) in the above three different limits, and by inserting the correspond- 



ing solutions in (IIII.32p . we study the behavior of S. 
I) In this case we can approximate A as 



X = ±l + -6\, |5A|<1, (III.33) 
6 



and equation flIII.30p as 

6X" + 36X' + 2(2 + 7)(5A ~ . (III.34) 

Then, solving the above equation, one obtains the following solution for A as a function 
of number of e-folds, at, 



SX ~ 6e-i°*(Ai cos(y ^ + 27 at) + A2 sin(y ^ + 27 at)). (III.35) 
As we see, the damping ratio is equal to — , ^ , which is less than one. Thus, SX 

' ^ ^ ^ 2^2(2+^' ' 

shows a damped oscillation and is damped within one or two number of e-folds. From 
the combination of (1III.35I1 and ( 1111.331) . we can determine A, which has the following 
form in the vicinity of A = ±1 



A ~ ±1 + e"i"*(Ai cos(y ^ + 27 at) + A2 sin(y J + 27 at)). (III.36) 

We see that A is exponentially approaching ±1 and the trajectory meets the attractors 
A = ±1. 

Inserting (1111.360 in (1111.320 . we can determine S 

which implies that S and 6X have opposite signs. 

As we see, for A^ in the vicinity of one a has the following behavior 

I — |<e-i°*, (III.38) 

here the subscript denotes an initial value. Thus is exponentially damped, with 
a time scale 2/(3d). 
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II) In the limit of very small values (|A| ^ 1) and considering the leading orders, 



equation flIII.30|) has the following form 

A" A' A'2 



— + 3— -3— -1-0, (111.39) 



which can be simplified as 



(^)" + 3(1)' + 21^0. (111.40) 



Solving the above equation, we obtain 



— ~ Aie-2<^* + ^26-^* , (III.41) 



which represents an exponential increase in |A| value with time scale of the order a~^. 
Thus, in the limit of initially very small A^ values, |A| is growing very rapidly and 
escaping quickly from the vicinity of zero. As a result, the above approximate solution 
is only applicable in first few e-folds where A^ is far from one. 

Eq. (IIII.4ip shows that for all possible values of Ai and A2, A monotonically increases, 
but interestingly this is not necessarily the case for a. In order to investigate this fact 
more precisely, we determine S for two different initial conditions in which (i) Ai = 
and (ii) A2 = 0. 



(i) Putting = in (IIII.4ip . we have A = Aoca"*, then after inserting in (IIII.32P 
we have 

S ~ -^!a2 -^e"*, (III.42) 
2 2' ^ ^ 

i.e. |S| exponentially increases in time. (Note that tjj = ^ is a constant in the 
leading order of e.) As mentioned before, due to the exponential growth of A, the 
approximation 0III.4ip is only applicable for the first couple of e-folds. Therefore, 
the growth of |S| can happen only during the first few e-folds, after that A gets 
close to one and |S| is exponentially damped. Note that although the phase of 
exponential growth of |S| lasts for the first couple of e-folds. 

(ii) Putting = 0, we obtain A = Aoe"*, leading to 



which is quickly damped. 
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Note that in this hmit ^ 1 and S has always a negative sign. 



Ill) In the limit of large values (|A| ^ 1) and recalling (IIII.23p . <^ 1, we have 



A^7 ^1. As a result, up to the leading orders, we obtain the following approximation 
for fUTLSOD 

A" + 3A' + 2A ~ 0, (III.44) 

which is identical to (1III.40I) that governs the evolution of p- in the limit of |A| ^ 1. 
Thus, the behavior of A in the limit of |A| ^ 1 is identical to the behavior of p- in the 
limit of |A| ^ 1, however, for completeness we present a more detailed analysis of this 
equation. 

Solving the above equation, we obtain the following solution for |A| ^ 1 

A ~ + Aae""*, (III.45) 

which has an exponentially damping behavior. Thus, in the limit of initially very 
large A^ values, |A| is damped very strongly and after one or two number of e- folds 
A becomes close to one and the approximate equation ( 1111.441) is not applicable any 
more. For all possible values of Ai and A2 in (IIII.45|) . A monotonically decreases, but 
interestingly this is not necessarily the case for &. To study this fact more explicitly 
we determine S for two initial conditions in which (i) Ai = and (ii) A2 = 0. 



(i) Putting y4i = in ( 1111.451) . we have A = Aqc which after inserting in (IIII.32I) 
yields monotonically increasing S 



Note that due to the exponential damping behavior of A, the approximation 
(1III.45|) is applicable only for at most one or two number of e-folds. Hence, the 



growth of S in this case, happens during the first few e-folds and after A gets 
close to one, E is exponentially damped. 

(ii) On the other hand, putting A2 = 0, we obtain A = Aqc"^"*, which gives S as 

S ^ ^ ^ !iAe-4"*, (IIL47) 
3 3 

which is quickly damped. 

In this limit A^ ^ 1 and the sign of S is always positive. 
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To summarize, assuming a system which undergoes quasi-de Sitter inflation in the sense 
that e is very small, we determined A and S. We see that regardless of the initial A values, 
all solutions converge to A^ = 1, within the few flrst e-folds. Note that, A^ = 1 corresponds 
to two values A = ±1, which are the isotropic solutions. As we saw before, A can not pass 
through zero during its evolution, so its sign does not change in time. As a result, as we have 
shown analytically and will be demonstrated numerically in the next subsection, system's 
trajectory eventually meets its attractor solution 

A ^ 1 if Ao > 0, 

A ^ -1 if An < 0. 



Furthermore, comparing (lIII.40p and (jIII.441) which are the approximate forms of (]III.30P in 



the limits of very small and very large A® values, we flnd out that the behavior of A in the 
limit of |A| ^ 1 is identical to the behavior of ^ in the limit of |A| ^ 1. 

Although A^ evolves toward the FLRW isotropic solutions, it is shown that in some 
solutions, grows rapidly at the first few e-folds saturating our upper bound of S ~ e for 
a short time. However, this growth stops fast (within a couple of e-folds) and S is damped 
for the rest of the quasi-de Sitter infiation. 

C. Numerical Analysis 



As seen from the action (jll.ip . our gauge-fiation model has two parameters, the gauge 
coupling g and the coefficient of the term k. On the other hand, the field degrees 
of freedom consist of two scalar fields ip and A, the isotropic expansion rate a and the 
anisotropic expansion rate a. Thus, our solutions are specified by eight initial values for 
these parameters and their time derivatives. The gravitational equations, however, provide 
some relations between these parameters. Altogether, each infiationary trajectory may be 
specified by the values of six parameters, {g, ipo, tpQ, Aq, Aq), here subscript indicates the 
initial value. 

In what follows we present the results of the numerical analysis of the equations of motion 
(1III.14I) . (1III.7I1 and the field equations corresponding to (p and A, for four sets of parameters 
corresponding to different positive initial A values (Aq = 1.05, Aq = 0.9, Aq = 0.1 and 
Aq = 10). Note that n, cto, ctq and ipo are given in units of Mpi. 
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FIG. 2. The classical trajectory for k = 1.733 x 10^^, g = 2.5 x 10~^ = 0.034, tpo = -10"^°, 
Ao = 1.05, Ao = 1.5 X 10^^*^. These values corresponds to a trajectory with do = 7.5 x 10~^, eo = 
3 X 10^3 and do = -8.6 x 10"^. 



► Discussion on diagrams in Fig^ and Fig[3i 

Figures [2] and [3] show the classical trajectories of two systems initially close to the isotropic 
solution (A = 1). In other words, these systems start from Ao values in the vicinity of A = 1, 
which one of them has Ao > 1 and the other one has Ao < 1. 

The top left figures in Figures [2] and [3] show classical trajectories of the field ip versus aot, 
while the top right figures show e{t). As we see there is a period of quasi-de Sitter inflation, 
where ip, e remain almost constant and e is very small. During this time, our initially 
anisotropic system mimics the behavior of the quasi-de Sitter inflation in the isotropic gauge- 
fiation [2|. Then, toward the end of the quasi-de Sitter inflation, e grows and becomes equal 
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to one (the top right figures), quasi-de Sitter infiation ends and ip suddenly falls off and starts 
oscillating. As we see from the top right figures, the parameter e has an upper limit which 
is equal to two. This is understandable recalling flIII.21|) and that p,^ is positive definite. 

The bottom left and right figures show evolutions of our two dimensionless variables A and 
S during the first few e-folds respectively. As we see in the left bottom figures, A follows an 
underdamped oscillation and quickly approaches one . Furthermore, from the right bottom 
figures we learn that is quickly damped. After the quick damping of the anisotropics within 
a few e-folds, system almost follows an isotropic gauge-fiation trajectory. 

► Discussion on diagrams in Fig\^ and Fig^ 
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FIG. 4. The classical trajectory for k = 3.77 x 10^^, g = W'^, = 0.6 x lO'^, i^o = 10"^°, 
Ao = 10, Aq = —3.6. Here A = — dA which corresponds to Ai = case in pil.44p and as ex- 
pected by our analytical calculations has a period in which S is positive and increasing. These 
values lead to a trajectory with do = 4 x 10~^, eg = 0.24 and do = —5 x 10^^. The initial value of 
e is rather large, but with in one number of e- folds it decreases and reaches 10^^. Note that value 
of e at the point of maximum S is equal to 0.05 (S ~ ^e), almost saturating our upper bound for 
anisotrpy S. 

Figures H] and show the classical trajectories of two systems initially far from the 
isotropic solution (A^ = 1). These systems start from two positive large and small A values, 
Ac = 10 and Aq = 0.1. 

Our analytical study reveals that although A evolves quickly toward the isotropic attrac- 
tors A = ±1, there are two independent solutions for the evolution of S in both of A^ ^ 1 
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FIG. 5. The classical trajectory for k = 1.99 x 10^°, g = 10"^ Vo = 0.099, tZ-o = 10"^°, 
Aq = 0.1, Aq = 1.288 X lO^'^. Here A = aA which corresponds to the ^2 = case in (|III.40|) . 
As has been predicted by the analytical calculations, S is negative and monotonically damped. 
These values give a trajectory with do = 4 x 10"'^, eg = 0.16 which is a rather large initial e value. 
As we learn from the top right figure, e decreases within a few number of e-folds to 0.5 x 10^^. For 
these values ctq = —0.8 x 10~^. 



and ^ 1 limits. In the first one, there is a short period of a few e-folds in which S rapidly 
increases and then decreases in time. On the other hand in the second case, S exponentially 
decreases. Due to the similarities between the two behaviors of both limits and in order to 
make the long story short, we chose to plot a system of the first kind for Aq ^ 1 limit, and 
a system of the second kind for Aq ^ 1 limit. In Fig. HJ a system with A = 10 and A = —aX 
is presented, which from the analytical results we expect to have a period of growing S. On 
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the other hand, Fig. [5] shows a system with A = 0.1 and A = dA, which we expect to have 
a monotonically decreasing S term. 

The top left figures in Figures H] and [5] show classical trajectories of the field ip with respect 
to aot, while the top right figures indicate dynamics of e. As we see there is a period of 
quasi-de Sitter infiation, where ip remains almost constant and e is almost constant and very 
small. During this time, our initially anisotropic system mimics the behavior of the quasi-de 
Sitter inflation in the isotropic gauge-flation [2|. Then, toward the end of the quasi-de Sitter 
inflation, e grows and becomes one (the top right flgures), quasi-de Sitter inflation ends and 
suddenly falls off and starts oscillating. As we see from the top right flgures the parameter 
e has an upper limit which is equal to two. This is understandable recalling ( ]III.2ip and 
that Pk is positive deflnite. 

The bottom left and right flgures show evolutions of our two dimensionless variables A 
and S during the flrst few e-folds respectively. 

In the left bottom flgure of Fig. HI we see that A which started from Aq = 10 quickly 
decreases and gets close to one. The right bottom flgure indicates that S, which is initially 
equal to 1.2 x 10^'^, shows a phase of rapid growth and saturates our upper bound S ~ e. 
More precisely, the peak value of S is Sl^^^^^. = 0.016 which is about |e|tp^^^. After its sharp 
peak, E decreases quickly and within few e-folds becomes negligible. At that point, system 
mimics the behavior of isotropic inflation. 

The left bottom flgure of Fig. [5] shows A, initially equal to Aq = 0.1, quickly evolves 
towards one. The right bottom flgure indicates S, which is initially equal to —2 x 10~^, is 
exponentially damped and becomes negligible. As a result, after a few e-folds, the system 
undergoes an essentially isotropic quasi-de Sitter inflation. 



IV. CONCLUSION 



We have studied nonlinear stability of the quasi-de Sitter gauge-flation with respect to the 
initial classical anisotropics. Assuming a system which undergoes quasi-de Sitter inflation, 
in the sense that e, the dimensionless isotropic expansion rate time variation (III.4p . is very 
small we examined the behavior of gauge-flation with respect to the anisotropics of the 
homogeneous background. We showed both analytically and numerically that this system 
has two attractor solutions A = (^) ^ = ±1 which regardless of the initial values of A, all the 
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solutions converge to them within a few e-folds. Here ipi and 1^2 are our field values. These 
two attractor branches, which are physically identical due to parity and charge conjugation 
invariance of our gauge-flation action (lll.ip . correspond to the isotropic quasi-de Sitter 
solutions. Thus, gauge-fiation is globally stable with respect to the initial anisotropics. We 
note that this is not a trivial result because in our system we have turned on vector gauge 
fields in the background and that the vector perturbations have a non-trivial source |2|. 

We showed that this attractor behavior happens for generic initial values of anisotropy 
A and the convergence to the isotropic attractor point happens very fast, within a few e- 
blds. This behavior may be contrasted with the anisotropic infiation model discussed in 
111 , [isl , where the model exhibits an anisotropic attractor with non-zero, but nevertheless 
small S ~ e at the end of infiation. 

Due to the fast damping of anisotropics in our model and that they do not last up to 
five to ten e-folds before the end of infiation, gauge-fiation predicts no detectable features 
of statistical anisotropy in the CMB temperature-temperature correlation function. In this 



respect gauge-fiation is simi 



the models discussed in 



11 



ar to the standard scalar-driven infiationary models and is unlike 



Hl3|- 



Furthermore, we found an upper bound on the value of |S| = ^ ~ e. Our numerical 
analysis reveals that in the extreme limits of ^ 1 and A^ ^ 1, there is the possibility 
to saturate our upper bound for a very short lapse of time. However, after reaching this 
maximum value, S is exponentially damped and soon becomes negligible. 

Here we focused mainly on the Bianchi type-I homogenous, anisotropic background met- 
ric. However, we expect our result, the global stability of the quasi-de Sitter gauge-fiation 
against the initial anisotropics and its isotropic attractor fixed point, extend over the other 
Bianchi type models. (Similar analysis for the case of power-law anisotropic infiation model 



15| has been carried out in 18|.) This expectation is based on the fact that the global 
stability observed here is a result of the pk dominance (c/. discussions in the opening of 
section [11]) and that this term can be written in terms of the single "effective inflaton" field 
(j) {cf. (111.141) . (111.101) ). These features still hold for general Bianchi backgrounds. We will 
expand on this argument in upcoming work fl9|. 

Our analytic and numerical computations shows that for very large and small initial values 
of A there is a region where anisotropy grows exponentially for a very short period, before 
getting exponentially damped to the isotropic fixed point. Although in these cases the system 
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does not follow the strict dynamics indicated by the Wald's cosmic no-hair theorem [g] for 
the very short period of time, the anisotropies are indeed damped within few Hubble times, 
in accord with the cosmic no-hair conjecture. In more general viewpoint, one can show that 
in general inflationary models do not strictly obey dominant energy condition assumption 
of cosmic no-hair theorem (c/. discussions in the second paragraph of the introduction) and 
hence in principle there is a possibility of violating this theorem. This possibility can be 
realized and can be more pronounced in the inflationary models involving vector (gauge) 



fields. In 



19| we will extend cosmic no-hair theorem such that it is also applicable to the 



cases like our gauge-flation and the anisotropic inflation model 
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